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Abstract. The gravitational waves generated in the coalescence of massive 
binary black holes will be measurable by LISA to enormous distances. Redshifts 
2: ^ 10 or larger (depending somewhat on the mass of the binary) can potentially 
be probed by such measurements, suggesting that binary coalescences can be made 
into cosmological tools. We discuss two particularly interesting types of probes. 
First, by combining gravitational- wave measurements with information about the 
universe's cosmography, we can study the evolution of black hole masses and 
merger rates as a function of redshift, providing information about the growth of 
structures at high redshift and possibly constraining hierarchical merger scenarios. 
Second, if it is possible to associate an "electromagnetic" counterpart with a 
coalescence, it may be possible to measure both redshift and luminosity distance 
to an event with less than ~ 1% error. Such a measurement would constitute 
an amazingly precise cosmological standard candle. Unfortunately, gravitational 
lensing uncertainties will reduce the quality of this candle significantly. Though 
not as amazing as might have been hoped, such a candle would nonetheless 
very usefully complement other distance-redshift probes, in particular providing 
a valuable check on systematic effects in such measurements. 



PACS numbers: 95.85.Sz, 98.80.-k, 98.80.Es 



As other contributions to these Proceedings will make clear, there are currently 
major uncertainties in our understanding of the astrophysics of massive binary black 
hole coalescences. We are quite unsure how often black holes merge, or what mass 
spectrum describes mergers, or how these quantities are likely to evolve as the universe 
evolves. What is certain is that if such mergers occur for binaries whose total masses M 
are roughly in the range 10^ Mq < {l + z)M < 10'' Mq, LISA will measure these waves 
out to redshifts of order 10. Even if these events are rare, LISA will measure them; and, 
because the source of these waves arises at such large distances, these measurements 
have the potential to provide detailed information about the large scale structure of 
the universe. In this contribution, we discuss the potential of such measurements as 
cosmological probes. 

Coalescing binaries can be considered standard candles because general relativity 
predicts a unique form for the two polarizations of the binary waveform. For LISA 
measurements, the strongest I — 2, m — 2 harmonic of the emitted waves has the form 

h+ - 1 + (L • cos m)] , 
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hx = -^-j—iTTfit)]'/' [l ■ h\ Sin Mt)] . (1) 

In this equation, M — rn'^^rri!^^ I[m\ + 7712)^^^ is the so-called "chirp mass". The 
vector L gives the orientation of the binary (it is the direction of the orbital angular 
momentum) ; h is the direction to the source according to an observer that rides along 
with the LISA antenna. The phase function <I>(t) depends on intrinsic parameters such 
as masses and spins, and so should be written <I>(f; mi, TO2, 5*1, 6*2). The frequency 
jit) = d^/dt, and r is the distance to the source. Note that the waveform Q does 
not include certain effects such as Lense-Thirring precession of the binary's orbital 
plane and multipoles other than I = m = 2\ this will likely impact the quantitative 
details of our results somewhat. 

By accurately measuring the evolution of the phase function we measure 

the binary's intrinsic parameters. In particular, the chirp mass M can be measured 
with exquisite precision: 5^A j M. ~ 10~^ is a reasonable expectation. This is because 
M. largely determines the total accumulated phase in the measurement, and hence 
is very sensitive to that number; cf. Refs. The LISA antenna pattern and 

modulations that arise from the antenna's orbital motion constrain L ■ n fairly well. 
The only remaining parameter in the waveform is the source distance r. Most error 
in determining r comes from correlations with orientation and position errors 4 ; in 
practice, r is likely to be measured to a precision Sr/r ~ 1 — 25% _5 . 

Interpreting these formulas becomes somewhat more complicated when the source 
generating and hx is at a distance where cosmological effects are important. 
Without going into the details (see 013 for further discussion), Eq. Q| still works 
provided we redshift all frequencies and timescales, and replace the naive distance 
measure r with the transverse comoving distance Dm (see Ref. |H] for a detailed 
description of cosmological distance measures). Putting these two replacements into 
Eq. and using the fact that the luminosity distance Dl = (1 + z)Dm, we find 

Ul L J 

hx - ^i^^i#^k/(t)]'/^ \l ■ n] sin m)] ■ (2) 

Now $(t) is strictly the measured gravitational- wave phase function (and f{t) is 
likewise the measured instantaneous frequency). It depends on redshifted values of 
the intrinsic parameters: 4>(t) = (1-|-2;)toi, (l-|-z)m2, (1 -|-z)^5i, (l-|-z)^52]. The 
reason for redshifting these parameters can be simply explained using dimensional 
analysis. In general relativity, a mass m can only impact the evolution of the system 
as a timescale = Gm/c^. (General relativity has no intrinsic scale, so the scales 
seen in any particular problem must follow from that problem's specific parameters, 
such as masses.) When the system is placed at redshift z, the timescale is redshifted. 
Thus, the apparent mass likewise picks up the factor 1 + z. Similarly, a spin 5' 
impacts the system as a squared timescale t| = GS/c'^ and picks up a factor (l-f z)^. 
As a consequence of this, the phase evolution of a cosmologically distant binary is 
indistinguishable from a "local" binary with redshifted masses and spins. 

Measuring the waves from a distant binary black hole coalescence thus provides us 
with two particularly interesting types of information — redshifted masses of the form 
(1 -|- z)m, and the luminosity distance (as well as information about the source's 
location on the sky and the spins of the binary's members). There are two direct ways 
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to exploit this information for cosmological studies. First, we can assume that we 
know the universe's cosmography. This allows us to build a map between z and Dj^. 
From our inferred z(Di,), we can break the mass-redshift degeneracy and learn about 
black hole masses as a function of redshift. Second, if it is somehow possible to obtain 
the redshift independently, one can use the simultaneous measurement of z and Dl 
to improve the cosmography. (In principle, there is also a third track: if one knows 
the mass spectrum of the binaries, or has enough events to statistically sample the 
range of the spectrum, then it should be possible to infer the cosmological properties 
of this distribution; see Ref. [H] for details. We thank Sam Finn for pointing this out 
to us. Given the large uncertainties in our understanding of the massive black hole 
merger rate and the likely breadth of the merger mass spectrum, we will not discuss 
this third direction here.) 

Let us begin by considering the first possibility. We will describe the universe 
using the currently popular "concordance cosmology" , with matter density given by 
flm — 0.35, dark energy in the form of a cosmological constant (that is, with equation 
of state p = wp and w — —1) with rt\ ~ 0.65, and a Hubble constant whose present 
value is Hq = 65 km/sec Mpc~^. We will assume the universe is precisely flat, and 
that the relative errors in I^a and Hq are ~ 10% ^Hl- These assumptions allow us to 
build a map from redshift to luminosity distance |Hj which is simple to invert (at least 
numerically). We then ask: How well can LISA measure the parameters characterizing 
massive binary black hole systems, particular the luminosity distance, redshift, and 
masses? To do this calculation, we use the restricted second post-Newtonian waveform 
described in Ref. and a description of the LISA antenna as described in Ref. 
see Ref. [H] for complete details. 

An example of a binary that is measured particularly well is shown in Fig. 
^ The histograms describe parameter measurement accuracies found by randomly 
distributing 100 binaries over the sky at z = 1 with random orientations, and with the 
merger time randomly distribed within an assumed 3 year LISA mission. Each binary 
is taken to have masses mi = TO2 — 10^ Mq. In many cases the redshifted masses 
and luminosity distance are measured with exquisite precision. The distribution for 
the distance peaks at 6Dl/Dl — 2%, and most events have SDl/Dl < 20%; the 
peak for the redshifted chirp mass is at 6A4z/-Mz — 10"'', with most events having 
SMz/A4z < 10~^; and the peak for the redshifted reduced mass is at Sfj.z/P'z — 4%, 
with most events at Sfiz/p-z < 10%. The redshift does not appear to be determined 
as well, but this is entirely due to our assumed errors in the cosmography. If the 
parameters determining the mapping between z and were known precisely, we 
would find Sz/z ~ SD^/Dl. Again, we emphasize that these distributions are 
computed with the restricted waveform given by Eq. JQ), are likely to change when 
effects we have neglected are taken into account. 

The example shown in Fig. is particularly good, but is not far off what 
can be achieved for a broad range of system masses in the rough band 10** Mq < 
(l-|-2:)Mtotai < (several) x 10® Mq. In this band, the redshifted masses and the distance 
are typically measured with a relative error of a few tens of percent. These numbers 
are also approximately independent of mass ratio, at least for mxjm^ > 1/10 or so 
— the loss in signal-to-noise ratio that comes from the reduced mass ratio is mostly 
compensated by an increase in the number of measured cycles, so that measurement 
precision remains roughly constant. 

We thus conclude that, at the very least and using gravitational- wave information 
alone, LISA will be able to provide useful and interesting data on black hole masses 
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Figure 1. Parameter measurement accuracy distributions for binaries witli nii = 
m2 = 10^ Mq at 2: = 1. In this case, the parameters are measured particularly 
well: the peak determination of luminosity distance is at SDj^/Dj^ ~ 2%, the 
peak in the redshifted chirp mass is at SA4z/Mz — 10~*, and the peak in the 
redshifted reduced mass is at S/i^/fiz ~ 4%. The relatively large error in redshift 
determination {5z/z ~ 15%) is because our cosmological model assumes that the 
cosmological parameters are themselves only accurate to about 10%. If those 
parameters were known precisely, we would find 5z/z ~ SDj^ /^L ■ 



and redshifts, making it possible to study the mass and merger history of black holes 
in the universe. Particularly at moderate to high redshift, this could provide a wealth 
of data on the formation and evolution of the universe's structures |51 1111 [T^ . 

Before moving on to the next track in our study — redshift determined 
independently — we would like to comment on the accuracy with which distances 
are determined. Distance error is strongly correlated with the errors with which the 
orientation and sky position of the source are determined: the measured waveform 
(which is a weighted sum of the two polarizations /i+ and hx) takes the form 

/imcas = [(l + ^)-^]'^" [^/]2/3^(angies) COS [$(t) + 0(angles)] . (3) 

The functions J-'(angles) and ^(angles) schematically indicate the dependence of the 
measured amplitude and phase on a source's position and orientation angles; see @] 
for further discussion and details. To measure the distance accurately, we must nail 
down these angles as precisely as possible. 

As has already been mentioned, the various source angles are determined 
by exploiting LISA's orbital motion. Roughly speaking, the "angles" indicated 
schematically in Eq. (|3J) are effectively time dependent from the viewpoint of an 
observer who rides along with the LISA antenna (though of course they are constant 
with respect to the solar system's barycentre). This time dependence modulates /imcas, 
with the exact modulation encoding the values of the source angles. As a rough rule 
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Figure 2. Comparison of distance determination for binaries at 2 = 1 with nii = 
m2 = 10® Mq. In the top panel, we assume that LISA'S noise gets extremely bad 
below 10"* Hz; the bottom panel assumes the cutoff is at 3 X 10~^ Hz. Distance 
determination improves by about a factor of ten in the lower panel. This is 
because the wider band allows LISA to follow the binaries' phase evolution for 
a much longer time: in the top panel, the binaries radiate in band for about 15 
days; in the bottom, the time in band is about 10 months. Controlling the low 
frequency performance will have an important impact on LISA'S science on high 
mass binaries. 



of thumb, measuring the angles weU requires that LISA move through at least one 
radian of its orbit, translating to a rough minimum of 2 months of observation to 
adequately pin down the source angles. Sources that don't radiate in band for long 
enough (typically high mass systems) tend to determine these angles badly and hence 
have poor distance determinations. A cure for these systems is to open the LISA band 
by reducing noise at the low frequency end. Figure|21compares how well the luminosity 
distance is determined under the assumption LISA's noise becomes extremely large 
below / = 10""* Hz (top panel) and below / = 3 x 10~^ Hz (bottom panel). Both cases 
look at measurements of binaries that have mi — m2 — 10^ M© and z = 1; as in Fig. 
^ we randomly distribute the binaries' orientations, sky positions, and merger times 
(within an assumed 3 year mission). In the first case (10^^ Hz cutoff), the binaries 
only radiate in band for about 15 days. Distance determination is concomitantly 
poor, with a peak in the distribution at SD^/Dl — 25%. (This actually isn't too 
bad, largely because these binaries generate a very strong signal.) When the cutoff is 
lowered to 3 x 10^^ Hz, the binaries radiate in band for 10 months, and the distance 
determination is dramatically improved — the peak is now at SDl / Dl — 2%. Keeping 
the low frequency behavior under control is clearly very desirable in order to study 
high mass binaries. 

Even in the best case, the position determination that LISA is likely to achieve 
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Figure 3. Comparison of distance determination for binaries at 2 = 1 witii 
mi = 1712 = 10® Mq. The top panel shows the distribution for determination with 
gravitational waves alone. The bottom panel assumes that an "electromagnetic" 
counterpart provides an independent pointing solution for the binary. This 
breaks many degeneracies in the parameter determination, improving the distance 
accuracy by about an order of magnitude (note the different scales in the two 
panels). 



is not great by the usual standards of astronomy — the best angular resolution 
is on the order of several to several tens of arcminutes. // it is possible to 
associate the gravitational waves from a binary black hole merger with some kind 
of "electromagnetic" counterpart, the situation improves dramatically. By getting an 
independent pointing solution for the binary, many of the degeneracies that impact 
distance determination are broken. We illustrate this in Fig. |21 The top panel shows 
the distribution of SDl/Dl that can be expected when only gravitational waves are 
used to analyze the binary. The lower panel shows how the distribution changes if 
a counterpart exists and provides an independent pointing solution. The distance 
precision is improved by about an order of magnitude in this case. 

Associating a counterpart with a merger is not going to be easy. The number of 
galaxies in each LISA "pixel" can reasonably be expected to number in the hundreds. 
Also, it is far from clear what kind of electromagnetic signature will characterize 
a counterpart. Some work discussing the kinds of counterparts one can imagine 
has appeared ; we hope that the promise of coordinating electromagnetic 

observations with gravitational-wave measurements will motivate additional work in 
this vein. 

Associating a counterpart with a gravitational-wave merger measurement offers 
another exciting possibility: from such an association, it may be possible to directly 
measure the event's redshift, rather than inferring it by combining the distance 
with cosmographic information In principle, this could provide simultaneous 
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Figure 4. Likelihood contours (1 sigma) for tlie matter density f!m and the 
equation of state parameter w (relating the pressure and density of dark energy, 
p = wp). We assume that the universe is fiat, and that the underlying cosmology 
has Vim = 0.3, ui = — 1. We compare parameter determination following 
measurement of two gravitational- wave candles (at 2 = 1 and 2 = 3), and 3000 
Type-la supernovae with the SNAP {21] (evenly distributed from z = 0.7 to 
z = 1.7). Gravitational lensing has a dramatic impact on the gravitational- wave 
candle, puffing the likelihood contour out by a large amount. The event rate of 
supernovae is high enough that lensing effects can be averaged out. 



measurements oi Dl and z, each with precision < 1%. Such a measurement, though 
possibly rare and difhcuh to make, would provide invaluable information about our 
cosmography (complementing other well developed probes such as Type- la supernovae 
|17l I18|). At the very least, because the systematics are so different from that of 
other candles, a merger would increase confidence in all candles (assuming that their 
measurements are in accord!). At best, because the merger candle could be exquisitely 
precise, it could impact cosmological parameter determination with high weight. 

In fact, as has recently been shown '19', all candles have a fundamental limit to 
their precision set by gravitational lensing, arising from the propagation of radiation 
through the lumpy, inhomogeneous universe in which we live. Lensing impacts 
gravitational waves exactly as it impacts electromagnetic radiation. A lens with 
magnification will cause an event whose true luminosity distance is to appear to 
be at distance Dhj ^JJi.- (Note that /i can be less than 1 — "demagnification" is in fact 
quite likely.) By convolving this error with the expected magnification distribution 
P(/i) Pni, we find that 5Dl/Dl ~ 5 — 10% due to lensing is probable (with some 
dependence on the redshift of the merger event) [20]. Although intrinsically of high 
quality, the actual effectiveness of the merger plus counterpart standard candle will 
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be significantly reduced by lensing. 

Figure^illustrates cosmological parameter determination using standard candles. 
We assume that the universe is flat, and that the total density is given by matter 
(flm) plus "dark energy" with equation of state p = wp. The figure shows the 1-cr 
likelihood contours in the Om-w plane for several cases. The heavy black line shows the 
contour that would be obtained if two merger events with counterparts are measured, 
neglecting gravitational lensing. This contour is nearly identical to (indeed, slightly 
tighter than) that expected for 3000 Type-la supernovae measured by the proposed 
SNAP satellite (dotted line) The dashed line illustrates what happens to the 

heavy black line when the systematic uncertainty induced by gravitational lensing is 
taken into account. It's rather sobering to note how much of an effect the lensing has 
on the gravitational-wave measurements. Lensing does not impact supernovae nearly 
as much: the supernova event rate is high enough that measurements can average 
away lensing effects, essentially sampling the full range of the lensing probability 
distribution. 

Combining these lensed gravitational wave events with the supernovae changes 
the supernova contour very little. We believe that, in the end, the most important 
contribution of a gravitational-wave standard candle will be as a check on systematic 
and evolutionary effects in the candle dataset. As a candle with drastically 
different properties, coordinated gravitational wave/electromagnetic measurements 
would improve confidence in all standard candles. 
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